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^ , We study the anomalous magnetic moment of the muon, a^, and lepton flavor violating decay 

O ' — 67 in TeV scale B — L extension of the Standard Model (SM) with inverse seesaw mechanism. 

■ We show that the B — L contributions to are severely constrained, therefore the SM contribution 

remains intact. We also emphasize that the current experimental limit of BR{^ 67) can be 
satisfied for a wide range of parameter space and it can be within the reach of MEG experiment. 



I. INTRODUCTION 



Ph. 
Oh' 

^ ■ The anomalous magnetic moment of the muon has been measured at Brookhaven National Laboratory to 

I , a precision of 0.54 parts per million. The current average of the experimental results is given by [ij 

in ! a^'^P = 11659208.0(6.3) x 10"^°, (1) 

\ which is different from the Standard Model (SM) prediction by 3.3cr to 3.6tT 0,0] 

2: Aa^ = a;^''P-a™ = (28.3±8.7to28.7±8.0) X 10-1°. (2) 

10 ! . . 

^ This discrepancy has been established by the impressive accuracy of recent theoretical and experimental 

results. Therefore, it is tempting to consider the above result as a strong signature for physics beyond the 

SM. It is important to note that the SM estimation for depends on the low-energy hadronic vacuum 

polarization, which is the main source of the uncertainty. The above result is obtained when hadronic 
- 1 — I 

I vacuum polarization is determined directly from the annihilation of e~ to hadrons. However, if hadronic 
H ' T decays are included, substantially larger value for a^'^ is derived that reduces the discrepancy to about 



2.4(7 only 3. 

In addition, non-vanishing neutrino masses confirmed by neutrino oscillation experiments are one of 
the firm observational evidences for an extension of the SM. The simplest way to account for small neutrino 
masses is to introduce right-handed neutrinos into the SM, which are Majorana-type particles with very 
heavy masses. In this case, type I seesaw mechanism can be implemented and an elegant explanation 
for light neutrinos is obtained. Recently, it has been shown that TeV scale right-handed neutrinos can be 
naturally implemented in B — L extension of the SM Q, where three SM singlet fermions arise naturally 
to cancel the U{1)b~l triangle anomaly. Also, the scale oi B — L symmetry breaking can be related to 
supersymmetry breaking scale Ql, therefore, the right-handed neutrino masses are naturally of order TeV 
scale. 

In order to fulfill the experimental measurements for the light neutrino masses with TeV scale right-handed 
neutrino, a very small Dirac neutrino Yukawa couplings, < 0(10^^) must be assumed Q. In this case, 
the mixing between light and heavy neutrinos are negligible, and hence the interactions of right-handed 
neutrinos with the SM particles are very suppressed. In Ref. ^ , a modification to the TcV scale B — L model 
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is proposed to prohibit type I seesaw and allow another scenario for generating the light neutrino masses, 
namely the inverse seesaw mechanism [lo| . In this scenario, the neutrino Yukawa coupling is no longer 
suppressed and can be of order one. Thus, the heavy neutrinos associated to this model are quite accessible 
and lead to interesting phenomenological implications. 

In this paper we analyze the anomalous magnetic moment of the muon in TeV scale B — L extension 
of the SM with inverse seesaw mechanism. We provide analytical formula for loop contributions due to 
the exchange of right-handed neutrinos, B — L gauge boson, and extra Higgs. We show that right-handed 
neutrinos give the dominant B — L contribution to a^. However, the unitarity violation limits of the light 
neutrino mixing matrix restrict this effect significantly. We also consider the impact of the right-handed 
neutrinos on the Lepton Flavor Violation (LFV) decays /i — 67. We show that the rate of this decay is 
enhanced and becomes within the reach of present experiments. 

The paper is organized as follows. In section 2 we briefly review the TeV scale gauged B — L model with 
inverse seesaw mechanism. We focus on the neutrino sector and show that the unitarity violation limits of 
Umns mixing matrix constrain the mixing between light and heavy neutrinos. In section 3 we study the 
anomalous magnetic moment of the muon due to the exchange of heavy neutrinos, B ^ L gauge boson Z' 
and B — L extra Higgs H' . In section 4 we analyze the LFV process /x — )- 67 and the constrained imposed 
by the experimental limit of BR{fi — e-f) on the heavy neutrino contributions. Section 5 is devoted for the 
numerical results and possible correlation between and BR{p — > 67). Finally we give our conclusions in 
section 6. 



II. TEV SCALE B- L WITH INVERSE SEESAW 

In this section we briefly review the TeV scale B — L extension of the SM with inverse seesaw mechanism, 
which has been recently proposed in Rcf. This model is based on the gauge group SU{3)c x SU{2)l x 
U{1)y X U{\)b-l, where the U{1)b-l is spontaneously broken by a SM singlet scalar x with B — L charge 
= — 1. Also a gauge boson Z' and three SM singlet fermions lyR- with B ~ L charge = — 1 are introduced 
for the consistency of the model. Finally, three SM singlet fermions 5*1 with B ~ L charge — —2 and three 
singlet fermions S2 with B — L charge = +2 are considered to implement the inverse seesaw mechanism. 
The B — L quantum numbers of fermions and Higgs bosons of this model are given in Table HI 
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TABLE 1: B — L quantum numbers of fermions and Higgs particles 
The relevant part of the Lagrangian in this model is given by 

+ {D^cf>)^{D^^) + {D^x)Hd^x) - Vi(l),x) - (XeUeR + KU^r + Xsi>RXS2 + h. 

- ^,s^ix''s,^^^s^,x'S2, 



(3) 



where F' ~ d^Z^ — d,yZ' is the field strength of the U{1)b-l- The general expression for the covariant 



derivative is defined as 



D^ = d,- igsT-Gl - ig-Wl - tg'YB, - ig"YB-LZ'^, (4) 
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where g" is the U{1)b~l gauge coupHng constant. The last two terms in Cb-l are non-renormahzable 
terms, which are ahowed by the symmetries and relevant for generating small mass for Si and 52 at TeV, 
are required by inverse seesaw mechanism. Few remarks are in order: i) The B — L symmetry allows a 
mixing kinetic term F^^F'^^^ . This term leads to a mixing between Z and Z' . However due to the stringent 
constraint from LEP 11 on Z — Z' mixing, one may neglect this term. In our analysis we assume a minimal 
model oi B ~ L extension of the SM. ii) In order to avoid other possible non-renormalizable term that may 
spoil the inverse seesaw mechanism that we adopt, a discrete symmetry like Z/i is imposed. Hi) In order to 
avoid a large mass term msSiS2 in the above Lagrangian, one assumes that the SM particles, vr, and 
^2 are even under a Z2-symmctry, while is an odd particle. Finally, x) is the most general Higgs 
potential invariant under these symmetries and it is given by Q 

x) = + mlx^X + M4>Uf + \2{x^xf + A3(x^x)(0V), (5) 

where A3 > — 2\/AiA2 and Ai, A2 > 0, so that the potential is bounded from below. 



The non- vanishing vacuum expectation value (vev) of = u'/\/2 is assumed to be of order TeV, 

which is consistent with the result of radiative B ~ L symmetry breaking found in gauged B — L model with 
supersymmetry (3]. The vev of the Higgs field 0: = v/\/2 breaks the electroweak (EW) symmetry, 

i.e., V = 246 GeV. After the B — L gauge symmetry breaking, the gauge field Z' acquires the following mass: 



The experimental search for Z' LEP II 
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M'z' - 9"'v'- (6) 
leads to 

Mz'/g" > 6 TeV. (7) 



Also, after the B ~ L and the EW symmetry breaking, the neutrino Yukawa interaction terms lead to the 
following mass terms: 

C'^ = mDVLVB. + MnVbS2 + h.c, (8) 

where itid = "^^^v and Mm = '^^sv' . In addition the second non-renormalizable term in Eq.Q induces 
a Majorana mass for S2 fermion. Hence, the Lagrangian of neutrino masses, in the fiavor basis, is given by 

/IJ'n = llsS2S2 + [niDVLVB + A/jvi^^5'2 + h.c), (9) 

where /is = -fjjj 10~^ GeV, hence M is of order intermediate scale 10^ GeV and the flavor indices are 
omitted for simplicity. Therefore, the neutrino mass matrix can be written as A4,^tp'^'ip with -ijj ~ (i'£, vb, S2) 
and Mu is given by 

' TOD 

M,^\ ml Mn \- (10) 

Mj; Us 

The diagonalization of this mass matrix leads to the following light and heavy neutrino masses, respectively: 

m^i = tudMn^ fJ-siMMy^ml, (11) 

™L = ^Ih' " + (12) 

It is now clear that the light neutrino masses can be of order eV, with a TeV scale if fig ^ M^. 
Therefore, the Yukawa coupling A^ is no longer suppressed and can be of order one. Such large coupling is a 
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crucial for testing this type of models and probing the heavy neutrino physics at LHC, as shown in Ref. 12|. 
The li ght neutrino mass matrix in Eq. (jlip must be diagonalized by the physical neutrino mixing matrix 



'MNS 



13| . i.e., 



Thus, one can easily show that the Dirac neutrino mass matrix can be defined as : 



mo = Umns V mtT'^ Ry l^s^ Mn, 



(13) 



(14) 



where R is an arbitrary orthogonal matrix. It is clear that this expression is a generalization to the expression 
of m£) in type I seesaw, which is given hy rriD = Umns \/ ^t]'^^ R \^Mn 



141. Accordingly, the matrix V that 



diagonalizes the 9x9 neutrino mass matrix Mu, i.e., V^M^V = A^^"*^, is given by [H 

V 

where the matrix V3X3 is given by 



^6x3 ^6x6 



V3, 



1 1 

2 



Umns ■ 



(15) 



(16) 



It is clear that in general Vsxs is not unitary matrix and the unitarity violation, i.e., the deviation from 
the standard Umns, is measured by the size of ^FF'^ . The matrix V3X6 is defined as 



V3 



3x6 



(03x3,i^)"^^6x6, F^moMJ,^ 



(17) 



Finally, Vgxe is the matrix that diagonalize the {vr, 6*2} mass matrix. Note that due to the Higgs mixing 
term in the potential V{(t), x)\^^ physical Higgs scalars {H, H') are given as a linear combination of and 
X, with the following masses 



(18) 



A detailed analysis for the phenomenology of the Higgs bosons of this model at the LHC has been considered 



Ref. 



If 



III. B-L CONTRIBUTIONS TO THE MUON ANOMALOUS MAGNETIC MOMENT 

In this section we analyze new contributions to the muon anomalous magnetic moment due to the extra 
particles of the B — L TeV scale model. From the effective Lagrangian of leptonic sector, one finds the 
following interactions 

c = ^{v;,D,rw+PLfi + v^.ni'^w-PLV^) + <?"/i7"^> + \ sine^H'^i, (19) 

where y is 9 x 9 extended MNS matrix, as discussed above, A^ is the Yukawa coupling of the muon, and 



d is the mixing angle between the SM-like Higgs and extra Higgs [16|. Thus, one can easily observe that 
the new contributions to the anomalous magnetic moment of the muon are generated by one loop diagrams 
involving the exchange of W gauge boson and heavy neutrino, or Z' boson and /i exchange, or H' neutral 
scalar boson and /i, as shown in Fig. [TJ Therefore, one can define a^~^ as 

(20) 
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FIG. 1: The new contributions of the muon anomalous magnetic moment m B — L extension of the SM. 



The calculation of the first diagram in Fig. [T] leads to 



<-%£^En:.^-:/(-.J, (21) 

where r,y^ = {m^./MwY ^-^d given by 

, 10 - 43r + 78r2 - 49r3 + 4r'^ + ISr^ ln(r) 
/W = ^(137)4 • (22) 

In this calculation, we assume that {m^/AlwY — 0- For r ~ one finds that /(O) = 10/3, while if r 1 
then f{r) 4/3. Thus, in the SM this contribution implies 

Kr = ^£|xy = 3-89xl0-^ (23) 

where the mixing matrix V is given by the unitary Umns mixing matrix, i.e., I^^m*!^ ^ 1- '^^^ model 

with TeV scale B L, the 9x9 mixing matrix V is unitary, however the 3x3 mixing matrix of light neutrino 
is no longer unitary. In our analysis, we constrain ourselves with the following non-unitary limits for light 
neutrino mixing matrix [3l 

0.994 ± 0.005 < 7.0 x 10"'"^ < 1.6 x lO^^ 
|7V7V^|w| < 7.0x10-5 0.995 ±0.005 < 1.0 x lO^^ |. (24) 
< 1.6x10-2 < 1.0x10-2 0.995 ±0.005 

In this case, one finds that 0.99 < I^m^P — 1- Hence the SM-likc contribution is slightly reduced to 

3.851 X 10-3. Since 4/3 < /(r) < 10/3, one can easily see that 

1 D . ^ 

> itM + l- 

i—l 

Thus, the ratio R'j^ = ajl/(ap^'^ lies within the tiny range: 

0.994 < i?;; < 1, (25) 

which means that within TeV scale B ~ L extension of the SM, the discrepancy between the theoretical 
prediction of the anomalous magnetic moment of the muon and its experimental measurement remains 2.4ct 
as in the SM and another source of new physics is required to account for this difference. 
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Next, we consider the contribution of Z' to a^. From the second diagram in Fig. [TJ one finds the foUowing 
result: 

where we assume that r„i^ = (mfj,IMz'Y — Oi g{rm^) is given by 

I \ Z\l-Z) r,„,.^0 1 



Hence one finds that 

< ''^^{0\<it^{m^) -2.34x10- (28) 

This contribution is quite smah and one can neglect the effect of the Z' diagram. Finally, we consider the 
diagram of extra Higgs. We find that the corresponding contribution to is given by 



A„ sin0P 2ml 



a,, 



327r2 m]j, 



-Hr'), (29) 



where = [m^/mH'Y — is assumed and A^^ = m^/v ~ 0(10 ■*). The loop function h{r'^^) is given by 

3(l-r;,J4 • (30) 

Hence one can estimate this contribution, for rriH' — 0(100) GeV and sin^? = 1/a/2, as 

<' ^ lo-^^ (31) 

which is also quite negligible. Therefore, one concludes that the B — L contributions to can not account 
for the reported discrepancy between the theoretical and experimental expectations. It is worth noting that 
B — L contribution to a^, obtained from the loop diagram mediated by heavy neutrinos can be of order 
the SM contribution and has a significant effect if the mixing between light and heavy neutrinos is sizable. 
However this mixing is strongly constrained by several leptonic processes [isf . 



IV. /x^ 67 IN TEV SCALE B- L WITH INVERSE SEESAW 

We now consider the LFV decay /i — ?> 67 in the TeV scale B — L model with inverse seesaw mechanism. 
Many experiments have been designed to search for LFV processes, in particular fi — > 67. The current upper 
limit is given by the MEG experiment [3| 

BR{p 67) < 2.4 X 10"^^ (32) 

New experiments are expected to improve this limit by three order of magnitudes. It is important to note 
that the SM result for the branching ratio of /i — > 67, with neutrino masses as in Eqs. (|42|43p . is given by 

BR{fi 6,7)™ ~ 10"^^ (33) 

Thus, the observation oi fj, ^ ej decay will be a clear signal for physics beyond the SM. We perform the 



7 




FIG. 2: 67 in B - L extension of tlie SM. 



calculation of /i — > 67 due to the exchange of light and heavy neutrinos and W gauge boson inside the loop, 
as shown in Fig. [21 The amplitude of /i — !■ £.7, in the limit of me — > 0, can be written as 



Let us define the coefficient of the amplitude A as 



9 



2—1 



Therefore, the decay rate is given by 



r(^^e7) = ^|ap. 



Using the dominant decay mode of r(/i — > ez^i^) ~ 7n|^G|,/(1927r^), the branching ratio is given by 



BR{p. 67) 



r(/i — > ei/j/) 



Stt m5G| 647r 



(34) 



(35) 



(36) 



(37) 



where a = e^/47r ~ 1/137. From the experiment upper bound in Eq.([32]), one finds the following constraint 
on the light and heavy neutrino mixing: 



< 0.000149. 



(38) 



In case of extremely heavy right-handed neutrinos, the lepton mixing matrix V^xa is almost unitary. 
Therefore, the contribution of light neutrinos, which corresponds to i = 1,2,3, is almost zero. In addition, 
the contribution of heavy neutrinos is quite suppressed due to the small mixing between light and heavy 
neutrinos in this case (V^i ~ Vei ~ mn/Miy ^ 0(10"^)). Hence the above constraint is satisfied and 
BR{^j. 67) < 10" ^2 is obtained. 

However, within TeV scale inverse seesaw the lepton mixing matrix is non-unitary. Also the mixing 
between heavy and light neutrinos are not small, since mo /Mm ^ 0(0. 1). Therefore, the bound in Ea. (|38)) 
can be written as 



^ „ 3 9 

i=l j=4 



< 0.000149, 



where rj ~ {mu^ . /MwY' and a-s defined in Eq.(|17p. is given by 

^<^Mj = [(0' moM^^) Fexe] i(2).,-3 = ( 0' Umns\J mt]''^R\[^\ Vg 



(39) 



(40) 



l(2)j-3 
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Thus, for rj ^ 1, i.e. f{rj) = 4/3, one finds 



< 0.0000636, 



(41) 



which implies that {FF'^)2i.i2 10 ^. This bound can be easily satisfied, due to the constraints imposed 
on the off-diagonal elements of the non-unitary U]\{ns mixing matrix. 



V. NUMERICAL RESULTS 

In our model oi B — L extension of the SM with inverse seesaw mechanism, the relevant input parameters 
involved in the computation of the anomalous magnetic moment of muon are the following: 

1. Three right-handed neutrino masses. 

2. Three mass parameters. 

3. Three angles of the orthogonal matrix R. 

In fact, the form of the Dirac neutrino mass matrix m^j, given in Ea.(|14p. guarantees that we obtain the 
correct light neutrino masses and mixing matrix Umns- 

The solar and atmospheric neutrino oscillation experiments provide the following results for the neutrino 
mass-squared differences with best-fit values within la errors j20| : 



Ami, = (7.64t°:l^) X 10"^ eV^ 



lAr 



'311 



(2.45 ±0.09) X lO"-' eV 



Therefore, one finds 



7.64 X 10 
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2.45 X 10- 



"1^1 1 1 



(42) 
(43) 

(44) 
(45) 



with arbitrary m^, . If <C 7.64 x 10 '^cV^, the ansatz of hierarchal the light neutrino masses is obtained. 
In this case one gets the following the light neutrino masses: 



< 10"^ eV, 



0.008 eV, 



0.05 eV. 



(46) 



In our analysis, we adopt these values for the light neutrino masses and also assume that the neutrino 
mixing matrix is given by Eq. (|24p . From Eq. (|40p . one notices that the mixing element V^i, which plays a 
crucial role in the result of a^, can be enhanced if: (i) fisi ^ . , (ii) the orthogonal matrix R is maximally 
mixing. For example, if = 2.7 x 10"^ GeV, Mn, = 900 GeV, Ma,, = 1500 GeV, Mn^, = 1900 GeV, and 
the other parameters are fixed as in Fig. |31 then one finds the following extended MNS mixing matrix: 



V 



/ 


0.806 


-0.591 


0.001 


-0.008 


-0.008 


-0.012 


-0.012 


0.0001 


-0.0001 






-0.418 


-0.569 


0.701 


0.020 


0.020 


-0.012 


-0.012 


0.068 


-0.068 






0.417 


0.570 


0.701 


-0.020 


-0.020 


0.012 


0.020 


0.068 


-0.068 















0.707 


-0.707 


-0.001 


0.001 





















0.001 


-0.001 


0.707 


-0.707 

































-0.707 


-0.707 






-0.032 


-0.026 





-0.707 


-0.707 





















0.029 





-0.001 


-0.001 


-0.707 


-0.707 










\ 








-0.136 














0.701 


-0.701 


) 



(47) 
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FIG. 3: BR{fi -5> 67) versus Sm = Mjv^ - Mn^, for A/jv^ = 100,500,900 GeV from up to down, respectively. The 
horizontal dashed line refers to the MEG experiment Hmit of BR{fj, — >■ 67). The other parameters are fixed as 
follows: Mjv3 = 2000 GeV, fis^ = 10"^° GeV, fis2 = 10"® GeV, = 2.62 x 10"^ GeV, m^,^^ = 10"" GeV, m„,^ = 
8.5 X 10"^^ GeV, m,,,^ = 5.05 x 10"" GeV. 

l.i — . — . — . — ^ — . — . — . — ^ — . — . — . — ^ — . — . — . — ^ — . — . — . — I 




0.995- 



2.x 10^* 4.x 10-* 6.X10-'* 8.X10-* 1.x 10"^ 
/^.„(GeV) 

FIG. 4: The ratio between the SM and B — L contributions to the anomalous magnetic moment of muon, iij^, as 
a function of ^^3, for Mn, = (120, 125, 300) GeV, (400, 465, 800) GeV and (900, 1500, 1900) GeV from up to down, 
respectively. The other parameters are fixed as in Fig. [3] 



From this example, one notices that the elements V^i, i = 4, .., 9 are of order 0(0.01) which induce the SM- 
like contribution 3.846 x 10~^ and total contribution of order 3.863 x 10"^. Thus, the ratio R"^ = aj^/(ap^^ 
is given by R""^ = 0.994. 

In the Fig. [3l we present the BR{^ — > ej) versus Sm = Mjy^ — Aljy^ for Aljy^ = 100, 500, 900 GeV from up 
to down, respectively. Here we assume that Mn^ = 2000 GcV, fisi = 10"^" GcV, fis2 = 10~^ GeV, fig^ = 
2.62 X 10-9 GeV, m^,^ = 10"" GeV, m^,^ = 8.5 x lO^^^ QeV, m^,^ = 5.05 x 10"" GeV. Note that the 
BR{^ — 7^ 67) is not sensitive to the value of A/jVa- For other LFV processes like r — >■ 67 and r — > /X7, the 
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present experimental limits of their branching ratios are given by 21 1 

BR{t ej) < 3.3 X IQ-^, BR{t ^7) < 4.4 x 10"**. 

One can easily show that in our models these experimental bounds can be translated into the following 
constraints: 



i=l 



< 7.5 X 10"^ 



i=l 



< 8.6 X 10 



-3 



This implies that (FF^)3i,i3 < 1.5 x 10"^ and (FF^ )23,32 ^ 1-7 X 10 ^. As mentioned above, F = 
itidM^^ ~ 0(0.1), therefore these constraints are naturally satisfied in our model, for mo ^ C'(100)GeV 
and Mn ~ 0{TeV) and no constraint will be imposed. 

In Fig. H we plot R"^ as a function of /i,,,, for Mn, = (120, 125, 300) GeV, (400, 465, 800) GcV and 
(900, 1500, 1900) GeV from up to down, respectively. The other parameters are fixed as in Fig. [31 In 
addition wc vary the angles of the orthogonal matrix R from to tt. As can be seen from this figure, there 
is no significant difference between the SM expectation and the total result of g — 2 in TeV scale B — L 
extension of the SM with inverse seesaw. 



VI. CONCLUSIONS 



In this paper we have computed the anomalous magnetic moment of the muon in TeV scale B — L 
extension of the SM with inverse seesaw mechanism. The one loop contributions due to the exchange of 
right-handed neutrinos, B — L gauge boson, and extra Higgs have been analyzed. Wc showed that right- 
handed neutrinos may give a significant contribution to a^, however it turns out that it is quite restricted 
by sever constraints from leptonic processes. Therefore, the SM contribution in i? — L extension of the SM 
with inverse seesaw mechanism remains intact. Thus the discrepancy between the theoretical prediction of 
and its experimental measurement requires a different source of new physics, like supersymmetric models 
with minimal fiavor violation, which usually respect the LFV constraints. We also studied the impact of the 
right-handed neutrinos on the LFV decay fi — > ej. Wc have shown that the rate of this decay is enhanced 
and is reachable by the MEG experiment. 
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